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1. Introduction 



It is believed that the hydrolysis of the adenosine triphosphate (ATP) is the universal mech- 
anism providing the energy for diverse biological processes such as photochemical reactions, 
cross-membrane ion transfer and signal transduction, muscle contraction, cellular mobility and 
transports to mention just a few pQ. All these processes involve protein molecules as mediator 
of the long distance energy transfer. However, it is yet not clearly understood how this energy 
can be transported along the polypeptide chain at long distances, without being dissipated or 
dispersed. 

An early explanation of this problem, based on quantum mechanical model, has been pro- 
posed in the mid-seventies by Davydov and collaborators [21 El H] . The essence of their theory 
is the assumption that the energy released by the ATP hydrolysis can be captured by the pro- 
tein molecules where it excites the high-frequency amide-I vibrational mode (C = O stretching 
quanta or vibrons) of a peptide group. Due to the dipole-dipole coupling this energy delo- 
calizes between the adjacent peptide groups giving rise to vibron band states. Such collective 
excitations may coherently propagate along the polypeptide chain participating in such a way 
in the energy transfer. Nevertheless, as it was shown by Davydov and collaborators [U El El H] , 
their lifetime is to short which rules out the long range energy transfer via band-like mech- 
anism. Instead, they suggested that the energy losses of the vibron may be prevented by its 
(self-)trapping in the potential well created by the induced local distortion of the molecular 
crystals. So created complex entity, i.e. the vibron surrounded with the local lattice distortion, 
may propagate in a soliton form along the chain with minimal energy losses preserving its shape 
and velocity for a long time. 

According to the general theory of self-trapping (ST) phenomena [6], the character of 
ST states is determined by the mutual ratio of the values of three basic energy parameters of 
the system. These parameters are the exciton (electron, vibron, etc.) bandwidth 2\J\ and char- 
acteristic phonon energy hue (determined by the phonon cut-off frequency uq) characterizing 
respectively, the time scale of motion in excitation and lattice subsystems. The third parame- 
ter is lattice deformation energy Eg which measures the strength of vibron-phonon interaction 
[21 El El [TO]- Two limiting cases could be considered as well understood: large- and 
small-polaron limit [51 El El 13 [10]. The first one is reached when the excitation bandwidth 
and lattice deformation energy both highly exceeds the characteristic phonon energy i.e. in the 
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adiabatic strong coupling regime. In that case phonons are slow with respect to the excitation 
dynamics and form, essentially classical, large radius quasi-static potential well in which that 
excitation (vibron) is trapped. Contrary, in the non-adiabatic limit, when the characteristic 
phonon energy is larger than the vibron bandwidth, the quantum nature of the phonons plays 
a crucial role, and small polaron (SP) formation takes place. They are characterized by a 
short-range lattice distortion that follows polaron motion instantaneously. 

The vibron bandwidth in hydrogen-bonded macromolecules such as polypeptide chains is 
about 0.001 eV, while the phonon cut-off frequencies are typically about 10 13 Hz [TTI [T2"| [T5]. 
As a consequence, basic assumptions for the soliton formation are not fulfilled [H]. For that 
reason, it has been suggested by some authors that vibron self-trapping in hydrogen-bonded 
macromolecular chains might result in formation of a small-polaron, rather than a soliton 
[IH IS [16l [TTI EH HH] • The applicability of such concept is restricted to the strong coupling 
and non-adiabatic limit which imposes certain doubts on the whole SP concept since only 
the non-adiabatic condition is satisfied, while the problem of the strength of vibron-phonon 
coupling in these structures is still under the debate. 

Usually, the studies devoted to vibron excitation problem in macromolecular chains involve 
the value of the strength of vibron-phonon coupling constant extracted from an indirect com- 
parison between experimental and theoretical results [IE]. According to these estimations, the 
value of this constant in hydrogen bonded macromolecular chains ranges between 30 pN and 
60 pN [16j [20]. However, there are such estimations of that parameter, performed by the ab 
initio calculations for the formamide dimmer, which suggests that its value is significantly lower 
|21[ |22| [23] . This fact became by the motivation for some authors to define criteria to discrim- 
inate between the weak and the strong vibron-phonon coupling limits [TTJ [21] . According to 
[24J, the strength of the vibron-phonon coupling in a lattice of if -bonded peptide units at 
room temperatures belong to strong coupling limit, but uncertainty of the value of the coupling 
constant may leads to completely different conclusions. So, this problem is still unresolved and 
requires a further research. Additionally, some recent numeric [25] and theoretical [26] studies 
of the ST phenomena in hydrogen-bonded macromolecular chains indicate that the character 
of ST states can be qualitatively very different depending on the values of system parameters, 
temperature and, therefore its proper theoretical description requires an approach that goes 
beyond the conventional strong-coupling SP theories. Especially useful in that sense is the 
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concept of partial dressing, which is applied in a quite wide part of system parameter space. 
It relies on modified Lang-Firsov (MLF) unitary transformation. Such an approach is in close 
correspondence with the supplementary variational treatment of the problem by means of the 
Toyozawa ansatz [27]. 

In this paper we study the character of vibron states in a-helicoidal macromolecular chains. 
We supposed that vibron interacts with acoustic and optical phonon modes. The problem of 
the vibron interaction with optical phonon modes has been intensively considered in many 
publications in the framework of the Holstein model [TTJ [28j [29] , especially after the indications 
that in ACN the vibron interaction with optical phonon modes plays a crucial role for the vibron 
small-polaron formation [501 EI] • On the other side, the problem of the vibron interaction with 
acoustical phonon modes usually has not been considered. This problem, in our opinion should 
not be disregarded since acoustical phonon modes have quite lower energies in comparison with 
the optical ones, and consequently, can be easily excited in the whole temperature range. 

Our study is performed within the D=l dimensional model resembling the real a-helix 
structure. We assumed that, due to the dipole-dipole interaction, the vibron excitation could 
be delocalized either from a some structural element to its first-neighbor ones (from n-th 
to (n ± l)-th structural element), or from a some structural element to its non-neighboring 
(from n-th torn/ (n, n± l)-th) structural elements. This light generalization of simple one- 
dimensional model enables treatment of the vibron derealization along the direction defined 
by the covalent bonds in a-helix (corresponding to the vibron derealization between the first 
neighboring structural elements), while the vibron derealization along the hydrogen bonds in 
a-helix macromolecule corresponds to vibron derealization between each third neighboring 
structural elements in the present model. We restrict our attention to the single-vibron case. 

Our particular interest is to establish conditions when vibron moves in a coherent (via 
band states) or incoherent (diffusive motion by means of the random hops between neighboring 
peptide groups) fashion along such structure. Special attention will be paid to the examination 
of the phenomenon known as small polaron crossover, which may occur for the certain values 
of system parameters and in the certain temperature range. This feature is connected with 
known ambiguities in SP parameters, effective mass, for example, arising as a consequence of 
the sudden jump in the magnitude of the degree of phonon dressing. 



4 



2. Model Hamiltonian 



The system under the consideration consists of single vibron excited on n-th structural 
element of the macromolecule, which physical properties are affected by the thermal vibrations 
of the macromolecular chain. We suppose that vibron excitation can move along the chain from 
n-th to the (n ± l)-th structural element, either from n-th to the (n ± 3)-th structural element. 
The first possibility corresponds to vibron derealization along the covalent bonds, while the 
second possibility corresponds to vibron derealization along the hydrogen bonds in 3D a-helix 
structure. Due to the fact that the energy of the dipole-dipole interaction between second- 
neighboring peptide groups is smaller than the energy of the dipole-dipole interaction between 
the first-neighboring ones [12] . the vibron derealization from n-th to (n ± 2)-th structural 
element is not considered. Corresponding Hamiltonian of the vibration excitation and phonon 
subsystem can be written in the following form [8]: 

H = A^4A n -^J g 4(A n+fl + A l - g )+E^ s J^ 

n n q 

+^=J2^F q e^ h AtA n (B q + Bl q ). (1) 

V -i» q n 

In the Hamiltonian above the quantity A is the vibron excitation energy, Al(A n ) describes 
the presence (absence) of the vibron quanta on the structural element, which is positioned on 
n-th lattice site, B^(B q ) creates (annihilates) phonon quanta, and u q is the phonon frequency. 
The inter-site overlap integral J g characterize the vibron transfer between macromolecular chain 
structural elements: g = 1 in the case of the vibron transfer along the covalent bonds, and g = 3 
in the case of the vibron transfer along the hydrogen bond. Derealization along the covalent 
and hydrogen bonds will be examined separately, and therefore the summation over the index 
g was not included in (Tj[|). This refers to two practically independent models describing the 
vibron derealization along particular bonds. In such a way, the mutual influence of the vibron 
derealization along covalent bonds on vibron derealization along hydrogen bonds (and vice 
versa) has been disregarded. Widely accepted data for values of the hopping constant J 3 is J 3 = 
Jhb = 7.8 cm -1 , while for the hopping constant between different spines of hydrogen-bonded 
peptide units is J\ = J cov = —12 A cm -1 [TTJ [T5J [32] • The parameter h is a distance between two 
neighboring structural elements along the axis of the macromolecular chain. Finally, F q = F* q 
is the vibron-phonon coupling parameter, which is a function of the coupling constants \ g . 
Coupling constant accounts for the modulation of the vibron frequency placed on the n-th 
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structural element due to the external motion of the (n + g)-th structural element. The value 
of the coupling constant that determines the vibron-phonon coupling in the direction defined 
by hydrogen bonds (xhb) m a-helix structure is known [161 120]. However, the influence of 
the motion of the (n ± l)-th and (n ± 2)-th structural elements on the amide-I vibration on 
the n-th structural element is still unknown. According some authors [32], numerical values 
of these parameters satisfy following conditions: Xhb > Xi » X2- In this paper we assume 
that xi — 0, Xi — and X3 — Xhb- Consequently, in the case of the vibron interacting with 
acoustic phonon modes, F q becomes as, F q = 2iXhb\J^jj^- srn (3g/i) [US EZ], while in the case 
of the vibron interacting with optical phonon modes this parameter becomes: F q = Xhb\J 2~it~T- 
Phonon dispersion law is determined by the geometric structure of the macromolecular chain, 
mass of the structural elements from which the chain is built, as well as the values of the 
elasticity coefficients that characterize the stiffness of the covalent and hydrogen bonds in the 
chain. Due to the fact that in the case of a-helix secondary structure three intermolecular 
interactions is the minimal number of bonds required to have a stabilized structure [33], we 
assumed that each structural element is connected with its first-neighbors, second-neighbor 
and third-neighbors. The elasticity coefficient K\ of the covalent bonds, that connects nearest- 
neighbor peptide units, varies between 45 N/m and 75 N/m [34J, while the elasticity coefficient 
K 3 of the hydrogen bonds ranges from 13 N/m to 20 N/m [HI [T81 [321 El]- In contrast to the 
analysis carried out in the [32J , we take into account non-zero value of the elasticity coefficient 
K 2 of the second-neighbor bonds which presence can significantly affect on the nature of the 
phonon spectra. We vary this value from interval K 2 G [mini^, max^], where mini^ << K$ 
and K\ > max^ >> K3. Consequently, dispersion law for acoustic phonon modes have 
following form: 

2 UC X . 2 (qh\ AK 2 . 2 4K 3 . 2 (3qh\ 

u n = —rv sin — H — — — sin (ah) H — — sin . 

q M \2 ) M yH ' M \ ^ ) 

In the case of optical phonon case we used the dispersioneless approximation: u q = uc, where 
phonon cut-off frequency is u>c = 2 J K ^+ AI ^2+9K 3 j^gj Here, M is the mass of the peptide unit. 
For the mass parameter we adopted the following numerical value: M — 5.7 • 10 -25 kg jllj . 

In order to examine under which conditions dressed vibron excitations represent dynam- 
ically stable eigenstates of the system, we cross to small-polaron picture using MLF unitary 
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transformation operator [26] 

U = exp {- J= £ E / g e-^ ft 4 A n (fi_ 3 - fit) | 

and we rewrite in terms of new operators a n = UA n U^ (at = UA n W) and 6 g = UB q W 
(b q = UB^W), representing the dressed vibrons and new phonons in the lattice with shifted 

equilibrium positions of molecular groups. 

p* 

Here f q = o^-j-, < 5 < 1 is variational parameter, which characterizes the degree in which 
the vibron distorts the lattice and the lattice feedback on the vibron, i.e. vibron dressing. 

On the next stage, we apply a simple mean-field approach, which enables the explicit 
accounting for the influence of thermal fluctuations on SP properties and, in particular, the 
degree of dressing. In the present case, the temperature enters in our model Hamiltonian im- 
plicitly through phonon fluctuations around its new equilibrium positions. Thus, the explicit 
temperature dependence may be introduced here by the appropriate averaging of the trans- 
formed Hamiltonian over the phonon subsystem. In particular, we define (in wave number 
representation) an effective, mean-field Hamiltonian H as follows: 

H = H + H int , (2) 

where H = H v + H ph , H v — (H — H ph ) ph and H int = H - H ph - (H - H ph ) ph . The symbol 
( } h denotes the average over the new-phonon ensemble, which is in thermal equilibrium state 
at temperature T. In such a way we obtain 

Ho = J2 E9 sp(k)a{a k + J2^ q blb q , (3) 

k <l 

H ™t = iEEW-^l^+'S-EEU^I^^) ( 4 ) 

VJ' 9 k fci k 2 

-e~ W ^S klM ) + e~^ h (T klM (-g) - e~ w ^5 klM ) } , 
where a k = -^J2n elknh(1 n, and operators T kl ^ 2 (±g) look as 

T klM {±g) = ^Ee^-^LA, ^ere 9 n = e v^^^>, (5) 

n 

while mean values of T kl ^ k2 (±g) over phonon subsystem calculated as follows 
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The variational energy of the small-polaron band state is given by the formula, 

E 9 SP {k) = A - 1 £ {(/, + f- q )F q - hu q \f q \ 2 } - 2J g e~ w ^ cos(gkh) (7) 

q 

with quantity 

W A T ) = ±= E \U\ 2 (2* 9 + 1)(1 - cos(gqh)) (8) 

q 

denoting the vibron-band narrowing factors, which characterizes the degree of the reduction 
of the corresponding overlap integrals or equivalently the enhancement of the polaron effective 
mass. For obvious reasons, they are sometimes called the " 'dressing fractions'" or the " 'dressing 
parameters'". At last, v q = l/(e hu}q / kB ' r _ i) denotes the phonon average number. 

Partially dressed vibron states represent dynamically stable eigenstates of the system pro- 
vided that variational parameter corresponds to the minimum of the system free energy. Thus, 
we search for its optimized value minimizing the system free energy, T = —k B T In Tr{e~ H / fcsT }. 
According to Bogoliubov variational theorem (which states that the function 

•F b = -A; B TlnTr{e-^ T } + (H mt ) Ho 

determines the upper bound of the system free energy: J 7 < J~b), the procedure of the mini- 
mization of the J 7 can be replaced by the procedure of the minimization of the Tb- In above 
expression, (H int )H denotes the mean value of the interaction polaron-phonon term, over non- 
interacting polaron-phonon system H . 

Using vector basis {\ip(k, {n q })) = a k \0) V ib ® Yi q — |0) p /i|, V/c, Vn^}, one can obtain follow- 
ing form of Tb'- 

T% = -k B T In E e~ E spW /kBT - k B T^ln 1 gT . (9) 

k q 1 c 

It is useful to have in mind that only first term of T B depends on variational parameter 5. 

3. Results and discussion 

In order to analyze the character of the dressed vibron states in macromolecular chains, 
we examine the degree of the vibron-polaron dressing (i.e. its effective mass) and the vibron- 
polaron ground state energy. Note, the ground-state vector of the effective Hamiltonian Ho 
reads \ipGs) — a \ |0)„j& ® T\ q |0) g ? where fig |0) g is the phonon vacuum vector and k = ir/3h 
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(g — 3), or ko — (g — 1). This case is also the most relevant for the spectroscopy since vibron 
bandwidths in these media are very narrow and its dispersion in all practical examination of 
optical spectra can be neglected [3 El I2S1 HSU EZ] . As a consequence the ground-state energy 
corresponds to the lowest level of SP energy, i.e., Eq S = {ipcs \Hq\ ipas) = E 9 SP {ko). 

Further calculations are done in terms of two independent dimensionless parameters: adia- 
batic parameter B = |^ determining the character of the lattice deformation engaged in the 
polaron formation and coupling constant S = (where E B = — J2g ^fff~ * s ^ ne lattice de- 
formation energy), which determines the polaron spatial size and the degree of dressing. This 
parameter, which has originally introduced by Holstein [HI |9] has special importance in the 
strong-coupling non-adiabatic regime (S » 1 and B « 1) where it characterizes the SP 
narrowing of quasi particle band which decrease exponentially as e~ s . 

The expressions for the dressing fractions and normalized ground state energy are calculated 
for vibron which delocalizes along the covalent bonds, and for the vibron which delocalizes along 
the hydrogen bonds. We take into account the vibron which interacts with acoustical and optical 
phonon modes separately. The variational parameter S is calculated by the minimization of the 
Bogoliubov upper bound of the system free energy: = and d d ^ 2 B > 0. 

A normalized ground state energy (normalized on the characteristic phonon energy and 
measured from the vibron excitation energy level: £q S = E( i s A ) takes the value, 



The first term in (flUj) corresponds to quasiparticle binding energy, while the second one 
corresponds to the width of the quasiparticle energy band. In the weak coupling limit or for 
intermediate coupling strength, there exists quasiparticle energy band of certain width. Con- 
sequently, in that case quasiparticle is delocalized (it can move along the chain on the coherent 
band-like fashion). On the other side, in the strong coupling limit (i.e. for the full dressed 
vibron case when 5 = 1) energy band is quite narrow, and consequently Sqs is proportional 
to quasiparticle binding energy (TOj [H] . In that case quasiparticle is almost localized on the 
particular macromolecular structural element (i.e. on the particular peptide unit), and it can 
move by random jumps between neighboring peptide units. Functional dependence of normal- 
ized ground state energy on coupling constant becomes £qs{S) — —S. These simple conclusions 
may be very useful in analyzing of the obtained results and explanation of the nature of the 
dressed vibron states. 




(10) 
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The expression for vibron dressing fraction, in the case of the vibron which delocalizes along 
the covalent bond reads 

when vibron interacts with acoustical phonon modes, while the quantity 



\V(r) = 5 2 £coth( — ) (12) 



is vibron dressing fraction when vibron interacts with optical phonon modes. In ( TTTT) integrals 
I w (t) = i Jq ^3 coth sin 2 3x(l - cosx)dx and I E b = ^ Jo sin ? 2 3x dx are obtained by the 
virtue of the standard rule, jjJ2q ■■■ ~ > ^ 1-1% ■■■dq- 

In the case when vibron delocalizes along the hydrogen bonds, vibron dressing fraction 
becomes 



W {t) = 5 2 S^^, (13) 

lEb 

when vibron interacts with acoustical phonon modes, and in the case when vibron interacts 
with optical phonon modes it takes the form 



W{r) = 5 2 S coth ^^j. (14) 

The integral in equation (fT3|) is Iw( T ) = ~ Jo ^ coth sin 2 3x(l — cos3x)<ia;. The system 
normalized temperature is r = ou x = — and x = qh. 

Obtained results are illustrated on figures Fig. 1-6. On figures Fig. 1-4, the dressing fractions 
and ground state energy are presented by the set of '"adiabatic"' curves W(S) and £gs(S) 
(each curve corresponds to a particular fixed value of the adiabatic parameter B) for various 
values of B, system temperature T, for both types of the bonds between peptide units, and for 
both cases of the vibron-phonon interaction (interaction with acoustical and optical phonon 
modes). The adiabatic curves £eg(S) obtained by the standard LF theory are presented by 
the thin curves, while the ones, obtained by the variational approach are represented by the 
thick curves. 
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Figure 1: Dependence of the vibron dressing fraction W (first row), and the vibron ground 
state energy £ GS (second row), on the coupling constant S, for various B and T. The case of 
the vibron which moves along the covalent bonds and interacts with acoustical phonon modes. 
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Figure 2: Dependence of the vibron dressing fraction W (first row), and the vibron ground 
state energy £ GS (second row), on the coupling constant S, for various B and T. The case of 
the vibron which moves along the covalent bonds and interacts with optical phonon modes. 
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Figure 3: Dependence of the vibron dressing fraction W (first row), and the vibron ground 
state energy £ GS (second row), on the coupling constant S, for various B and T. The case of 
the vibron which moves along the hydrogen bonds and interact with acoustical phonon modes. 



13 




Figure 4: Dependence of the vibron dressing fraction W (first row), and the vibron ground 
state energy £ GS (second row), on the coupling constant S, for various B and T. The case of 
the vibron which moves along the hydrogen bonds and interacts with optical phonon modes. 



One should be noted that from the above presented results in the strong coupling limit and 
for all values of B and T the obtained values of both quantities are qualitatively very similar 
with ones derived by the usage of standard SP approach (i.e. with results obtained with using 
of the Lang-Firsov unitary transformation). This similarity increases with the decreasing of B. 
For those values of the system parameters both approaches predict heavy dressed vibron states. 
In addition, with the growth of the system temperature and for fixed value of B, the values of 
the parameter S for which our approach and standard SP theory yield similar results moves 
towards lower values. These observations seem to be not surprising, because of the framework of 
the standard small-polaron approach, the growth of the system temperature has practically the 
same effect as the growth of the coupling constant. This point appears by the consequence of 
the fact that the SP dressing fraction (and consequently quasiparticle effective mass) increases 
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as temperature grows, but it also increases with the growth of coupling constant value [38J . 

On the other hand, for small and medium values of S our results are significantly differs from 
ones predicted by standard LF theory. In contrast to the standard SP theory, in that limit our 
approach predicts slightly dressed or partially dressed vibron states (instead of characteristic 
parabolic shape, the adiabatic curves £qs(S) has linear dependence of coupling constant). For 
fixed value of system temperature and S this difference increase with increasing of B. Similarly, 
for fixed S and B the difference increases with increasing of T. All these conclusions are fulfilled 
for both types of the intermolecular forces between peptide units, and for both cases of the 
vibron-phonon interactions. 

As it is demonstrated on Figs. 1-4 there are two clearly distinguished areas in the system 
parameter space, where W(S) and £qs(S) display qualitatively different behaviors. In fact, 
for each T, type of the vibron-phonon interaction, and bond type there exists the critical 
adiabatic parameter Be, such that, for all B < Be, dressing fraction W(S) and ground state 
energy Sgs{S) changes smoothly as a function of coupling constant S. In that case function 
W(S) monotonically increases, while £cs(S) for small values of S slightly decreases, and after 
some value of S it becomes very closely to linear decreasing function. Thus, for B < Be 
the vibron states for small values of S corresponds to small-polaron band states with gradual 
transition towards the self-trapped ones as the case the coupling constant increases. 

But when B reaches some critical value Be, adiabatic curves W(S) have an infinity deriva- 
tive at a point with abscissa Sc- In the same time, ground state energy £es{S) ceases to be 
smooth at the point with the same abscissa. The point in the system parameter space, with 
coordinates (Sc,Bc) determines the "'critical point"'. By further increasing of the value of B 
(B > Bq, i.e. for supercritical values of B), adiabatic curves are smooth up to a supercritical 
value of S, where they suffer the abrupt change. The main characteristic of the system in this 
region of parameter space is the spontaneous abrupt transition between the low dressed and 
high dressed states, which takes place through the discontinuous change in the magnitude of 
SP parameters, effective mass and ground state energy, in particular. As one can remark, with 
increasing of adiabatic parameter B, the Sc moves to larger values. But, as system temperature 
T increases, the Sc moves to lower values (for fixed value of B). One can also remark that 
this observations are stronger expressed for vibron which delocalizes along the hydrogen bonds, 
especially in the case of the vibron interaction with acoustic phonon modes. 
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The existence of the critical and supercritical points is a consequence of the appearance 
of multiple local extremes of the function Fb{§)- In fact, in the case of supercritical values 
of B, and for small values of S, the function J-b{$) is increasing with respect to 5. With 
increasing of S it firstly appears one local minimum at a point which abscissa is 5 ~ 0. With 
further increasing of S, there appear two additional local extremes: one local maximum and 
one local minimum (new minimum appears at the point with abscissa 5 ~ 1). For smaller 
values of S, the minimum that corresponds to the lower value of 5 is deeper, and consequently, 
for these values of 5" and B vibron is lightly dressed (small value of 5 implies a small value of 
the dressing fraction). When S reach Sc, both local minima have equal depths. This means 
that there exist two vibron states (light dressed vibron state and heavy dressed vibron state) 
which both correspond to the same minimal value of Tb- With further increasing of S the 
deeper minimum of Tb becomes the minimum that corresponds to heavy dressed vibron state 
(5 ~ 1). This model predicts a similar vibron nature for both types of intermolecular chains, 
and for both types of vibron-phonon interactions. 

Above mentioned discontinuities are well known and sometimes have been considered as 
unphysical and the artifact of the variational methods mostly used in particular calculations. 
Nevertheless, in spite of all doubts the rapid increase of the polaron effective mass with the 
coupling constant in the adiabatic regime (B 3> 1) is the main characteristic of the so called SP 
crossover and represents the general conclusion regardless of the particular theoretical method 
which is used pH E7J |39j EUJ HU H21 03] . Qualitatively precisely the same behavior has been 
undoubtedly confirmed recently by means of the numerically exact diagonalization [H] of the 
Holstein model where discontinuous transition! from free electron to immobile small polaron 
m e ff = oo has been predicted. Moreover, the possible coexistence of free and localized states 
have been recently pointed out by Hamm and Tsironis [25J who had investigated the dimen- 
sionality effects on large to small polaron crossover in Holstein's model. The emergence of 
energy barrier separating the two types of solutions has been predicted irrespectively of the 
system dimensionality. This barrier is responsible for the double-peak structure of the wave 
function and its appearance has been interpreted as a clear signature of the coexistence of free 
(large) and localized (small) polaron solutions in a certain parameter range, which is in quali- 
tative agreement with our results. Finally, it is interesting that there exists the close analogy 



The notion of phase transition has been used in [33] to emphasize its abrupt character. 
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between abrupt change from delocalized to localized vibron state in macromolecular chains and 
thermally activated positronium state transitions in alkali halide cristals. It is known that the 
ground state of positronium in such structures is delocalized (Bloch-type), but there exists addi- 
tional metastable state which is localized and it can be achieved by the increasing of the system 
temperature. The localization of the positronium was investigated experimentally [I5 l H6l Wf\ . 
and theoretically H9] , and it can be explained by the thermally activated transition from 
ground state to excited metastable state. 

As we mentined above, the critical values of coupling constant Sc and adiabatic parameter 
Bq depends of the nature of the intermolecular bonds, phonon nature, as well as of the system 
temperature. This dependence is represented on Figs. 5 and 6 by dashed lines. Supercritical 
values of S and B are presented for a several values of T by sets of "isothermal" curves (full 
lines) . 




Figure 5: The dependence of the boundary values of the system parameters (S, B) on temper- 
ature in the case of vibron derealization along the covalent bonds: a) vibron interacting with 
acoustical phonon modes; b) vibron interacts with optical phonon modes. 
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Figure 6: The dependence of the boundary values of the system parameters (S, B) on temper- 
ature in the case of vibron derealization along the hydrogen bonds: a) vibron interacting with 
acoustical phonon modes; b) vibron interacts with optical phonon modes. 

Dashed line divides the system parameter space on two areas: for the system parameter 
values that belong to the left area the curves W(S) and Sgs(S) are smooth, and for larger 
values of S they are very close to the ones predicted by the standard LF approach. Here, 
presented model predicts a significant difference with standard LF approach only for medium 
values of S, where obtained adiabatic curves have a characteristic shape. The points that lie 
at right side from the dashed line are supercritical points, where W(S) and Sqs(S) have an 
abrupt change. 

The critical values of system parameters are shifted to lower values, with increasing of 
the system temperature. As an example, in the case of the vibron which delocalizes along the 
covalent bond and interacts with acoustical phonon modes, these values belong to the adiabatic 
and strong coupling region. They enters in non-adiabatic region approximately at T = 145 
K (at room temperatures T = 300i^, we obtain Be = 0.92 and Sq = 0.303. These values of 
the system parameters belong to non-adiabatic region of intermediate interaction. In the case 
when vibron interacts with optical phonon modes, Bq belongs to adiabatic region for whole 
temperature range (at T = 300 K, we obtain Be = 1.55 and Sc = 0.819). However, in the case 
of the vibron derealization along the hydrogen bonds, and interaction with acoustical phonon 
modes, Be belongs to non-adiabatic region for all values of the system temperatures (at T = 1 
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K we obtain Be = 0.82 and Sc = 1.179), and at room temperatures it enters deeply in non- 
adiabatic region (at T = 300 K we obtain B c = 0.39 and Sc = 0.056)! In the case of vibron 
interaction with optical phonon modes obtained results are very similar to ones, obtained for 
vibron derealization along the covalent bonds and its interaction with optical phonon modes. 

4. Conclusions 

Summarizing, the presented model of the a-helicoidal macromolecular chain predicts that 
vibron self-trapping would results in the creation of partially dressed small-polaron band states. 
With the growth of the temperature, heavy dressed, practically localized small-polaron states 
may appear instead of the band states, irrespective on the nature of the bonds between two 
structural units. The values of the critical parameters for which the vibron nature changes are 
dependent on the system temperature and nature of the bonds. It is noticeable that values of 
the Be decrease with the rise of the system temperature. 

In the case of an a-helix structure our system model predict the following. When vibron 
interacts with optical phonon modes, the critical value of adiabatic parameter belongs to adi- 
abatic region for all values of system temperature, as well as for both intermolecular bond 
types. In that case the critical value of coupling constant Sc belongs to strong coupling limit 
at low temperatures, and reaches to intermediate coupling strength approximately at room 
temperatures. 

When vibron interacts with acoustical phonon modes, Be enters into non-adiabatic region 
at T = 145 K in the case of vibron derealization along the covalent bonds. However, in the 
case of vibron derealization along the hydrogen bonds, Be belongs to non-adiabatic region for 
all values of system temperatures. 

Due to the fact that the values of the system parameters for the alpha-helix macromolecules 
belongs to the non-adiabatic and weak to intermediate coupling limit, it seems that at room 
temperatures in such systems may appear only heavy dressed vibron which moves along the 
covalent bonds, but for vibron which moves along the hydrogen bonds may appears weakly 
dressed vibron states. 

Obtained results are similar with results obtained in the study of the electron polaron 
effective mass around the V point of the Brillouin zone (121 S3] . Additionally, according to the 
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study by Hamm and Tsironis [25], it follows that the results obtained here are consistent, at 
least on a qualitative level, with the exact numeric ones. 

In the latter respect, our study represents the background for investigations of the polaron 
properties in the three-dimensional structure of a-helicoidal proteins as well as the bound- 
ary between coherent and incoherent polaron motion. According to some authors (32J [33] the 
phonon spectra in a-helicoidal structures exhibits two acoustical phonon branches and one op- 
tical phonon branch. The first ones in long- wavelength limit are essentially polarized along the 
axis of the direction defined by the axis of the helix (longitudinal vibration mode) and torsional 
mode. The second branch refers to a radial motion of the helix structural elements (breathing 
mode). Due to the fact that the vibron-phonon coupling parameter for the three-dimensional 
problem depends of the projection of the phonon polarization vector on the direction of the 
vibron motion, it is necessary to investigate the vibron interaction with all three phonon modes 
in an appropriate manner. In addition, we stress the further studies require the inclusion of 
the mutual influence of both dressing fractions and evaluation of the temperature dependent 
mean lifetime of the localized states. 
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